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We prove a generalization of the van Leeuwen theorem towards quantum electrodynamics, pro- 
viding the formal foundations of a time-dependent Kohn-Sham construction for coupled quantized 
matter and electromagnetic fields. We circumvent the symmetry-causality problems associated with 
the action-functional approach to Kohn-Sham systems. We show that the effective external four- 
potential and four-current of the Kohn-Sham system are uniquely defined and that the effective 
four-current takes a very simple form. Further we rederive the Runge-Gross theorem for quantum 
electrodynamics. 
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I. INTRODUCTION 

Time-dependent density functional theory (TDDFT) 
Q is a formally exact reformulation of nonrelativistic 
quantum mechanics in terms of the one-particle density 
instead of the wave- function. It is the time-dependent ex- 
tension of the highly successful density functional theory 
(DFT) 0. The Runge-Gross theorem [3[ proves under 
some restrictions, that all physical observables are func- 
tionals of the one-particle density alone. Hence, instead 
of the complex wave-function on configuration space one 
only needs the simple one-particle density to fully de- 
scribe a multi-particle system. Further, the van Leeuwen 
uniqueness theorem Q provides a way to calculate the 
one-particle density of an interacting quantum system in 
terms of a noninteracting auxiliary system of fictitious 
particles moving in a local effective external potential. 
This so called Kohn-Sham construction makes ab initio 
calculations for complex time-dependent many-body sys- 
tems feasible This is due to the fact, that the time to 
numerically solve a noninteracting Schrodinger-like equa- 
tion only scales linearly with the number of particles, 
whereas in the case of an interacting Schrddinger equa- 
tion, it scales exponentially Q. 

TDDFT has been extended, e.g., to open quantum sys- 
tems Q or to include special relativity. In the time- 
independent theory a generalization to relativistic sys- 
tems is well-known [1, 0| • A relativistic description of a 
quantum system becomes important, e.g., if we consider 
condensed matter systems with a high atomic number. 
The relativistic generalization of DFT towards the time- 
dependent domain is until now mainly used to calculate 
linear response properties of heavy elements [e.g. see Q 
and references therein]. The main problem here lies in 
the approximations of the functionals, which in practice 
always have to be made. The formal foundation of the 
generalization of TDDFT to relativistic systems is given 
in There, the most fundamental theory we have in 
describing interacting electrons, i.e., quantum electrody- 
namics (QED), is used to deduce a fully relativistic func- 
tional theory. A functional description based on QED 



should in principle be able to describe all possible effects 
we observe in condensed matter systems. Hence, any ex- 
tension of the usual density functional approach, e.g., to 
superconducting systems [lOj, should in principle be de- 
ducable from this theory. In 0] a Runge-Gross theorem 
for QED is formally proven. To fully describe the coupled 
matter and radiation quantum fields one is in need of the 
four-current and four-potential of the interacting systems 
as those are the conjugate functional variables with re- 
spect to the external time-dependent four-potential and 
the external time-dependent four-current defining the 
QED Hamiltonian. Therefore, any observable of the sys- 
tem can in principle be calculated by only knowing these 
conjugate variables. However, the action functional ap- 
proach to a Kohn-Sham-like system and the associated 
effective four-potentials as well as four-currents exhibit 
the same drawbacks as the corresponding nonrelativis- 
tic approach: it leads to a symmetry-causality paradox 

mm 

In Sec. II, we introduce the QED Hamiltonian and de- 
scribe the formal problem. Then, in Sec. Ill, we em- 
phasize the symmetry-causality problem of the action- 
functional approach to a Kohn-Sham-like description Q 
of the interacting system. In Sec. IV we will deduce 
the central theorem of this work, formally proving the 
uniqueness of a Kohn-Sham-like construction for QED 
systems. Further we can restate the relativistic Runge- 
Gross theorem of Q and show that the effective four- 
current of the Kohn-Sham system takes a simple form. 
Further we present the defining equation for the effective 
four-potential and give a straightforward approximation. 
Finally we conclude in Sec. V. 



II. QUANTUM ELECTRO DYNAMICAL 
DESCRIPTION 

We start our considerations by choosing the initial 
configuration of the combined matter-radiation system 
Po = Y,iPi\*i(to))(%(to)\, with Pl > and ^Pi = 1 - 
The density matrix operator po describes the initial con- 
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figuration of the massive particle system as well as the 
photon field on the corresponding Fock space. The dy- 
namics of the system is governed by the standard QED 
Hamiltonian in fl3l [T3 | , which in the Heisenberg picture 
(we indicate the explicit time-dependence of the external 
Hamiltonian by t) and in units such that h = c = 1 reads 

H(t) = Hm + He + Hi n t + -ffext(i)- (1) 
Here the normal ordered (::) Dirac Hamiltonian of the 
matter field in terms of the four-spinor operators ip( x ) = 
'ijy (x)7° , ip(x) and the usual gamma matrices is defined 
as 



Hyi = 



ip(x) ( — i 7 ■ V + mo ) tp(x) 



(2) 



The components of the spinors obey the equal-time anti- 
commutation relations 

{$ a (x, t) , $a % t) } = 7 V 3 (x-y)- (3) 
In contrast to Q the electromagnetic part in terms of the 

four-potential operators (x) and (x) is given by the 
Fcrmi-Hamiltonian UM 



d 3 x : l-A» (x)A li {x) + {d k A v (x) ) (8 k A v (x) ) :(4) 



with k = 1,2,3. We use quantization in Lorentz gauge 
[HQ, i.e., dnrlpoA^x)} =: d' 1 {A ll {x)) = 0, and hence, 
the four-potential operators are subject to the covariant 
equal time commutation relations 



An{x,t),A v %t) =ig liV S 3 (x-y) 



(5) 



with the standard metric tensor g^ u having the sig- 
nature (+,—,—,—). Note that one has to use the 
Fermi-Hamiltonian and hence the Fermi-Lagrangian C = 
— l/2{d il A 1/ {x)){d v A ll {x)) in order to be compatible with 
the above equal time commutation relations. However, 
one can show that the expectation value corresponding 
to the Hamiltonian He can be rewritten into the usual 
form LUI 



(H E ) 



d 3 a; (E 2 (x) + B 2 (x) 



(6) 



with 



E(x) = —(A(x)) — (WAo(x)), 
B(x) = (V x X{x)). 
The interaction Hamiltonian is 



Hint := / A A x e f{x)A^{x) , 
where the four-current density is defined by 

j**(x) = : ip(x)j ,1 ^(x) : 



(7) 



ip(x),j^ip(x) 



(8) 



and e is the electron charge. The last equality is shown 
to be true in [T3|. This immediately highlights the con- 
nection between the normal ordering and symmetrizing. 
Both make the current behave correctly under charge 
conjugation [13 , [l5| . Formally we can instead of the 
normal ordered Hamiltonian use the symmetrized ver- 
sion thereof [HI, [l6| ■ This form has certain advantages 
when renormalizing the Hamiltonian [16| . The deduced 
equations of motion will not depend on the form chosen. 
Finally the external Hamiltonian is defined in accordance 
to as 



#ext(*) = J d'xe [j^(x)a^(x)+j^(x)A ll {x)) . (9) 

All classical external four-potentials a® xt (a;) are assumed 
to obey the Lorentz gauge condition, i.e., <9 M a° xt (a;) = 0. 
Because we want the initial configuration to be fixed, 
the possible gauge-transformations of the four-potentials 
have to leave the initial configuration unchanged 
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Further, the external four-current densities obey the con- 
tinuity equation, i.e., d^j^ix) = 0. The external four- 
current is either used formally in field theory to gener- 
ate the different interaction contributions or in radiation 
source problems. In a Kohn-Sham approach to fully rel- 
ativistic systems, the effective external four-current will 
be important to generate the interacting four-current and 
four-potential in the noninteracting system. 

We are interested in solutions of the fully interacting 
field theoretical problem. However, in order to find fi- 
nite answers we need to renormalize the so-called bare 
quantities of the theory [lg. Otherwise the local inter- 
actions introduce divergences in the theory. These so- 
called ultraviolet divergences are usually eliminated by 
the renormalization program of QED [here we do not 
consider the infrared divergences as they can be taken 
care of, e.g. by enclosing the system in a large box with 
periodic boundary conditions [Ha]]. The first step in the 
standard renormalization procedures is to regularize the 
theory, i.e. to introduce some kind of high energy cut-off, 
or equivalently, a smallest length scale. This makes the 
theory finite, however, also dependent on this cut-off. In 
QED this is routinely done via dimensional rcgularization 
[14| preserving gauge invariance. Via perturbation theory 
one then redefines the coupling constant, the mass pa- 
rameter and the field operators order by order such that 
the divergences are cancelled by so-called counter terms 
and the dependence on the cut-off is removed [la ]. The 
perturbative renormalization to all orders for QED with 
a time-independent external potential has been closely 
examined in [l5j . where a Hohenberg-Kohn theorem for 
QED is formally deduced. For a time-dependent external 
four-potential the renormalization of QED is analyzed in 
[l8| . Assuming a stable vacuum, i.e. the external poten- 
tial is sufficiently weak not to close the gap between the 
positive and negative eigenvalues [l6|, one can deduce 
general renormalization rules. It is important to note 
that the renormalization factors are defined by vacuum 
QED, i.e. they do not depend on the external fields. This 
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fact permits the comparison of two systems with differ- 
ent external four-potentials. One can also reformulate 
the external-potential problem in QED in terms of the 
well-established Furry picture [l6j . Again the prerequi- 
site of a stable vacuum is crucial [l9| . We therefore 
only consider external electromagnetic fields well below 
the so-called Schwingcr limit or critical field of QED of 
E CT = 10 16 V/cm where the QED vacuum becomes un- 
stable to e + e~ pair-production [2(J. Having this con- 
dition we can (at least perturbatively) renormalize our 
problem at hand. This condition is in agreement with 
the caveat raised in Yet, also for the renormalized 
theory we do not know if a solution exists. This, how- 
ever, will be tacitly assumed. 

As in the derivation of the relativistic Hohcnbcrg-Kohn 
theorem [l5| one might wonder why we actually use the 
problematic field theoretical frame work in the following 
and do not use an approximate relativistic many-body 
scheme, e.g. Dirac-Coulomb-Breit approximation. As 
stated in [15j on page 366, in such an approach the 
renormalization is done via the no-pair approximation 
which spoils the gauge invariance of the theory. In order 
to avoid this we retain the field theoretical treatment of 
the problem. 

The renormalized Hamiltonian of eqn. (TTJ) is de- 
fined uniquely by the choice of the classical external 
four-potential (up to a gauge transformation) and 
classical external four-current, i.e., 



H(t)=H([a°?,j£ ct ]-,t). 



(10) 



Hence for every different set of these functional vari- 
ables, we find different solutions pit) = /5([a® xt , jg Xt ]; t) 
to the interacting problem after we have fixed the ini- 
tial configuration pq. In an analogous manner, the solu- 
tions |$(f)) of a Schrodingcr equation with only scalar 
external potentials v(r,t) are uniquely defined by the 
external potential up to a gauge transformation, i.e., 
|$([u];i))($([u];i)| = p([v];t) after we have fixed the ini- 
tial state |$o)- The density matrix operator-form cancels 
the phase ambiguity of the gauge. In the usual TDDFT 
it is shown, that the conjugate functional variable to v, 
that is the one-particle density, will describe the sys- 
tem uniquely. Here, the conjugate functional variables 
to (a® xt ,j'c X t) are, as can be seen from eqn. ©, the ex- 
pectation values of the four-current j M (x) := (j M (a;)) and 
the four-potential A^(x) := (A^(x)). Our goal is to show 
that the density matrix operator pit) is uniquely defined 
by the four-current and the four-potential of the coupled 
matter-radiation system, i.e., 



(ii) 



and hence (j^jA^) defines the system uniquely. This is 
the Runge-Gross theorem for QED [|[ . For a fixed initial 
configuration there is a bijective mapping (j 71 , A„) i— >• p(t) 
such that all observables can be calculated in terms of 



(j M , Ay), i.e., for a self-adjoint operator O 

(d(t)) = o([po,r,A^t) 



(12) 



is a functional of [j^,A^) and the initial configuration 
Pq. This leads to the possibility of calculating all observ- 
ables in terms of (j^, An) alone. On the other hand, we 
want to show that for every interacting problem associ- 
ated with the Hamiltonian (fTJ there is one and only one 
nonintcracting Hamiltonian, i.e., Hi nt = 0, with an ex- 
ternal four-potential a° ff and an external four-current jg ff 
coupling either to the matter system or to the photon sys- 
tem respectively, which leads to the same pair (j' 1 ,A A1 ) 
as the interacting problem. This means, that one can 
solve a quantum system of nonintcracting fermions sub- 
ject to the so-called effective external potential a® to- 
gether with a photonic system under the influence of an 
effective external source instead of the original in- 
teracting problem. This amounts to a Kohn-Sham con- 
struction for QED. However, as we will discuss in more 
detail at the end of this work, one can in general not sep- 
arate the fermionic part of the quantum system from the 
photonic one in the nonintcracting auxiliary system. If 
it is possible, then we have two separate equations, one 
for the four-current and one for the four-potential of the 
quantum system. Both separate solutions then consti- 
tute the solution of the fully interacting theory. 
Now we will describe the shortcomings of the initial at- 
tempt to formally define a Kohn-Sham construction in 
QED. 



III. SYMMETRY-CAUSALITY PARADOX 

In reference 0] an approach similar to Q is pursued, 
where an action functional of the form 



W,A M ] = B\j»,A„] 



(13) 



Co 



dt / d 3 x (r(x)a^(x)+j^ t (x)A^x)) 



with 



B\j^A, t }= / dt (V(t)\id t -Hu-Hv- J?i„t|tf (*)) (14) 

Jto 

is used. Here we have lost explicit covariance and hence 
have fixed a particular reference frame. In principle we 
could keep explicit covariance by using the Tomonaga- 
Schwinger approach [l6[ and general space-like surfaces 
a(x). Then the integration becomes 



t' p fa(x) 

dt / d 3 x -> / d A x 

to J J cr (x) 



(15) 



and the time-derivative becomes a functional derivative 

W)M ' |. w> ^ l»W>-l«W> ,(n) 

Oa(X) f2(aj)->0 U(x) 
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where f2(x) is the four- volume between the space- like sur- 
faces a and a'. Consequently, the Hamiltonian is re- 
placed by its density. However, such general space-like 
surfaces are not more physical than the one we have cho- 
sen initially, i.e., the one with t = constant. Therefore we 
give up the explicit covariance for the sake of a simplified 
notation. 

Note further, that for this variational approach, and 
hence also in this section, only pure states, i.e., wave- 
functions \^(t)), not general density matrix operators 
are considered. Therefore our approach, that will be pre- 
sented afterwards, is also a generalization with respect to 
the states possible. 

Following ^ the action-functional A should be station- 
ary at the exact time-dependent currents and fields, i.e., 



= and 



SA^x) 



= 0. 



(17) 



From this stationarity principle the corresponding Kohn- 
Sham equations are derived. However, this approach will 
lead to a violation of causality. To show this we will fol- 
low the same reasoning as in [l2j . 

Assume the variational principle (|17l) to be valid. Then 
we find by B = A + J (j"a« 



t+jeVW, where / 



8j»{x) 



a™* (2) and 



SA^x) 



Now we have a one-to-one correspondence between 
(j^, Ap) and (a^Si^t) due to the Runge-Gross theorem 
for QED systems 0. Hence we can use the Legcndrc 
transformation to derive 



Ao-T'&t] = ~B+ (i% Xt + Je M xt^) • (19) 



This leads to 



= j^ix) and 



A^x). (20) 



We then can define the linear response kernels of the 
matter and radiation systems as 



6 2 A[a™\ 



Jext] 



Sa'^ ct (x)Sa° xt (y) 
6 2 A[a™\& t ] 



WW 
Sa cxt (y)' 

SA^(x) 



(21) 



^ext( a; ) <5 iext(y) ^cxt(y) 

Those kernels describe how the system reacts in first or- 
der at space-time position x to an external perturbation 
applied at space-time position y. It is obvious, that these 
kernels should be causal, as the system should not react 
to perturbations in the future. However, the second func- 
tional derivatives in eqn. (|21l) are symmetric with respect 
to their arguments. Hence we meet a contradiction if we 
assume the variational principle of eqn. (|17[) to be true. 



A thorough analysis of this contradiction [l2j shows, that 
in contrast to variations in the wave-functions, for which 
the action-functional of eqn. (|T3l) is stationary at the 
solution point, variations in the external four-potential 
or four-current fix the variations of the wave-function for 
all later times. Therefore one is usually not allowed to 
pose the boundary condition \5^>(t')) = for the varia- 
tion of the wave-function at the endpoint of the variation 
interval t = t'. Further, also one can not vary the imag- 
inary and real part of the wave-function independently, 
as those are also uniquely defined by the variation of the 
external four-potential or four-current. Solutions for this 
causality problem in a variational approach to TDDFT, 
e.g., b y using the Keldysh time-contour, are found in 

Gf&fii. 

Now we will derive a van Leeuwen theorem for quantum 
electrodynamical systems. 



IV. VAN LEEUWEN THEOREM FOR 
QUANTUM ELECTRODYNAMICAL SYSTEMS 

In a first step we consider the dynamical behavior of 
the four-current density operator j M (x). The Hciscnberg 
equation of motion reads 



id t j»(x)= j"(x),H(t) 



(22) 



As 7 71 is not explicitly time-dependent in the Schrodinger 
picture the usual total derivative of the Heisenberg equa- 
tion of motion is equal to the partial time-derivative. Af- 
ter some calculations we arrive at: 



w t r(x) = m ^{x) [yy - 7 v] fa) 



(23) 



+ip{x) 



7 M 7° 



17 • V ) + ( -rf ■ V ) 7 C V 



4>(x) 



+efax) [ 7 M 7V ~ YW] ${x)(A„{x) + cf^ix) 

Further we are interested in the equation of motion of 
the four-potential operators, i.e., 

idtA^x) = [i„(x),#(t)] . (24) 
After some calculations we arrive at 

id t A tt (x)=il fl (x). (25) 
The second derivatives with respect to time are found as 
(idt) 2 A^x) = -d k d k A^(x) - ej^x) - ef xt (x), (26) 

which are just the usual inhomogeneous Maxwell equa- 
tions in Lorentz gauge, where k = 1,2,3. 
To deduce a relativistic van Leeuwen theorem we now 
assume the external potential a cxt {x) and the external 
four-current j^y±{x) to be analytic in time, i.e., 
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has a nonzero convergence radius, and analogously for the 
four-current. Further, we assume that the quantum elec- 
trodynamically calculated four-current j^(x) = (j M (x)} 
and the four-potential A^(x) = (A^(x)) are analytic in 
time as well. We define in accordance to the nonrelativis- 
tic van Leeuwen proof Q : 



+ mo|(!t) [7^7° ~ 7°7 M ] $(a:) (28) 



ft"" (a;) = e^(x) [yW - 7 l/ 7°7 M ] $(<b) 
With this, eqn. becomes 



(29) 



(30) 



(31) 



The next derivative with respect to time becomes due to 
the Heisenberg equation of motion 

(idtfr(x) = H(t)] + [CtO^ #(*)] (32) 

n^(x),H(t)} al xt (x) + (id t a° xt (x)) n^x). 
This can be rewritten as 

[\d t ff{x) = (id t q^(x)) + (id t C t (x)) (33) 
+ (id t n^(x)) al xt (x) + {id t al xt (x)) h» v {x). 

For all higher orders we find accordingly 



(id t ) l+1 j»(x) = (id t f[C n (x) + qUx)} 
Wf) mraT(x)) ({id t ) l - n h^(x)) 



(34) 



n=0 

If we define 



j3 M (x) = d k d k A fl (x), 



(35) 



an analogous equation holds for the time-derivatives of 
the four-potential operators 

(id t ) l+ % l (x) = - (id t ) L [Mx) + ej^x) + ej^(x)\ .(36) 

These are the defining equations for all orders of 
the Taylor-expansion of the four-currents and four- 
potentials. Hence we can construct the analytic 
expectation value of the four-current and four-potential 
by taking all orders of the time-derivatives at t = to and 
averaging with the initial configuration po- This can be 
done in a successive manner. 

For a noninteracting system with Hamiltonian H'(t) 
analogous to (pj except for 



(37) 



we find the corresponding defining equations 
(id t ) l+1 j'»(x) = (id t ) l q£ a (x) 



(38) 



+ 



1 

E 



((id t ) n af(x)){(id t y- n h'^(x) 



and 



{id t ) l+ %(x) = - (id t ) l [p'^x) +ejf(x)] 



(39) 



Now we want an analytic effective four-potential a° ff (a;) 
and four-current j^g{x) as well as an initial configuration 
p' such that for the expectation values of both systems 



j'»{x)=f{x) and A'(x)=A ll {x). 



(40) 



The initial configuration p' of the noninteracting system 
must fulfill 



f{xM) = tr pQ'Ufi 



(41) 



Afj,(x, t ) 



tr 



PoA'^x) 



Further, the noninteracting initial configuration has to 
fulfill 



A^{x,t ) = tr p' A'Jx) 



(42) 



Now, in order that both analytic currents are the same, 
all orders of their Taylor-expansions have to match. 
Therefore, the expectation values at t = to for eqn. (|3~4")) 
with po and for eqn. (j3"5)l with p' have to be the same. 
For I = we find a defining equation for a^(x,to) as 

n'^(x, t )af(x, to) = [q^ n (x, t ) + q? nt (x, t )-q^ lIL (x, to)] 



+n^(x,t )a B * t (x,to). 



(43) 



All we need to assume is that vl^ix, t ) ^ 0. This seems 
well justified in general. For the analytic four-potentials 
to agree we need that the expectation values at t = to for 
eqn. (|36p with po and for eqn. (|39p with p' coincide. For 
I = we find therefore a defining equation for j£ s (x,to) 
as 



ejeg(x, h) = [ej^ix, t ) + ej£ Kt *o)] 
+ \p^(x,to)- P " 1 (S,t )] 



(44) 



However, we have p ,l {x,to) = d k dkA^ (x , t ) = p'^{x,to) 
and therefore 



^s(x,to) = ^(x,t )+j^(x,to)} 



(45) 



We have now defined the zeroth order of the Taylor- 
expansions for the effective four-potential and the ef- 
fective four-current. Note, however, that the zero- 
component of the four-potential, i.e., ag ff (a?, to), is not 
uniquely defined via the above equation as n /#1I/ = for 
p, = or v = 0. The chosen initial configurations and 
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a choice of the zero-component of the zeroth-order Tay- 
lor expansion fix a certain gauge for the effective exter- 
nal four-potential. The higher order zero-components are 
successively defined via the Lorentz gauge condition. 
In a next step we require that (id t ) 2 j^{x)\ t= t = 
(idt) 2 j' ,1 (x)\t=t - Hence, we find accordingly 



n"*»(x) d t af (x)\ t=t = d t + g^(x) 

\t=t a +n^{x)d t aT{x)\ t=to 
+< xt (x) d t n^(x)\ t=tn - af(x) d t n'» v (x)\ t= 



(46) 



t=to ■ 

The time-derivatives of, e.g., q^ n (x), at t = to 
are found with help of the Heisenberg equation, i.e., 
{[q^i n (x), H(t)])\ t= t - Only the initial configurations and 
the afore-defined zeroth orders of the effective potential 
and current is needed. Analogously one finds from the 
expectation values of the four-potential operators for the 
interacting as well as the noninteracting system 



dtj^{x)\ 



(47) 



Hence, for all higher orders / we find, if we assume both 
systems to lead to the same current-density as well as 
four-potential expectation values 

n'^(x) d\af{x)\= d\ [q^(x) + q? nt (x) - q£ a (x)] | t=to 



l 

E 

n=0 
l-l 

E 

n=0 



t=t 



and 



°tJ e S 



( x )\ 



t=t 



t=t 



(48) 



The time-derivatives on the right hand sides include 
terms of the effective potential and effective current up 
to order I — 1. Hence, we can successively construct all 
orders of the Taylor-expansion and find the effective four- 
potential and the effective four-current for a noninteract- 
ing system which leads to the same expectation values 
of the current and potential operators as the correspond- 
ing interacting system. It is important to note, that the 
effective external current is just the sum of the exter- 
nal current and the current of the matter system, i.e., 
3cs( x ) = Joxt( x ) + 3^{ x )- There is of course the possibil- 
ity that the convergence radius of one of the constructed 
Taylor-series is 0. This would mean that the terms of 
the series increase faster than l\a l for a an arbitrary pos- 
itive constant [T3|- This seems rather unlikely and we 
will assume the constructed series convergent. Via ana- 
lytic continuation [3| we can then formally extend both 
Taylor-series to all times. With this we can formulate 
the relativistic van Leeuwen theorem. 

Theorem Let p be the initial configuration of an in- 
teracting matter- electrodynamical system. Assume, that 



for the analytic external four-potential a° xt (x) and four- 
current (x) the corresponding calculated four- current 
j IJ '{x) and four-potential A^(x) are analytic in time as 
well. 

If we find an initial configuration of a noninteracting aux- 
iliary system p' , subject to the constraints 



f(x,t ) = tr p' Q j'»(x) 



A^ix, t ) 



tr 



PoA'JS) 



tr 



we can formally construct an unique effective external 
potential a° ff (x) and an unique effective external current 

3 



if^{x), depending only upon the two initial configurations 



and the chosen external potential and currents, such that 

f{x)=j'»{x), 
A,(x) = 4 Or). 

The corresponding effective and external potentials are 
uniquely defined up to a gauge transformation. If we as- 
sume the second system to be interacting as well and start 
from the same initial configuration we find, that there is 
one and only one pair (a® xt ,j° xt ) up to a gauge transfor- 
mation leading to (J 1 *, A^). The gauge of A^ is fixed by 
the initial configuration chosen. So the gauge freedom is 
only with respect to the external potential and not with 
respect to the quantum system. This corresponds to the 
quantum electrodynamical Runge-Gross theorem . To- 
gether with the above derived van Leeuwen uniqueness 
theorem for QED systems, this leads to the possibility 
of calculating all interacting observables with the help 
of auxiliary noninteracting systems. Such a procedure is 
a quantum electrodynamical generalization of the time- 
dependent Kohn-Sham construction. Note, that in Q 
the effective four-current had an extra term (5B)/8A V , 
This term is zero, as shown in our derivation. 

It is important to stress not only the dependence of 
the effective external four-current and the effective ex- 
ternal four-potential but also the Kohn-Sham procedure 
as a whole on the initial configurations po and p' . As- 
sume the noninteracting initial configuration takes the 
form p' = p' M (to) <8> Pe(^o)- Then the expectation val- 
ues of the noninteracting system decouple into a mat- 
ter part and a photonic part. Further, the matter part 
Pm(*o) = X)ii>»|A(*o))(A(*o)| is assumed to only con- 
sist of single Slater determinants \Di(to)) of orthonormal 
orbital spinors ipj(x). Therefore we can treat this nonin- 
teracting problem via solving a set of Dirac-Kohn-Sham 
equations 0,0 

\d t ip 3 {x)=^a-[-\V -a e s{x)} + m a /3 + a° s (x)^pj(x),(49) 
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as the nonintcracting Dirac equations decouple into single 
particle equations, and a corresponding Maxwell equa- 
tion with the effective four-current [|| 

d u F^ = d^A v {x)-d v A^{x)) (50) 
= e.?e Xt (.T) +ef(x). 

Here, the Dirac-Kohn-Sham orbitals tpj(x) are used to 

construct the exact current j^ix) = Y] ■ Cjipj(x)a IJ 'ipj(x), 
where Cj is the occupation number of the j-th spinor 
orbital. Although usually one is only interested in the 
matter part of the quantum system, even in this spe- 
cial case an explicit calculation of the four-potential is of 
interest to deduce approximations for the effective exter- 
nal four-potential. It seems reasonable to assume that an 
approximation based on the potential mediating the in- 
teraction leads to good results. At the end of this chapter 
we will deduce such an approximation. Note further, that 
a similar set of equations, i.e., one for the matter system 
and one for the electromagnetic field, is in principle also 
needed, e.g., in the calculation of the response of solid 
state systems via TDDFT [23| . If we do not assume this 
special form for the initial configuration, then the Kohn- 
Sham procedure will be more complicated, and we can 
in general no longer use simple Dirac-Kohn-Sham equa- 
tions. Further, in general the expectation values have to 
be calculated with the full configuration of the matter- 
photon system and will not separate as in the case above. 
From the considerations leading to the van Leeuwen the- 
orem we also have a defining equation for the effective 
four-potential. Due to the assumption that the interact- 
ing and the noninteracting systems should lead to the 
same four-current we have with eqn. (|31[) for both sys- 
tems 

n'^(x)af(x) = q£ in ( x )+€t(x)-tfjx) (51) 

A straightforward approximation for the effective four- 
potential can be deduced as in the nonrelativistic case 
[24 |. if one assumes p(t) ~ p'{t) = p' M (t) ® p' E (t). Then 

af(x)^a^(x)+A,(x) (52) 

as the radiation field and the matter fields decouple and 
1hit( x ) = n ' f1 '' (x)A„(x). Note that this approximation 



is closely related to the self-consistent field approach for 
external-field QED [Hj]. 



V. CONCLUSION 

In conclusion, we have formulated a quantum elec- 
trodynamical van Leeuwen theorem circumventing the 
symmetry-causality problem of the action-functional ap- 
proach to relativistic Kohn-Sham systems. Instead of 
solving the interacting matter-photon problem, one can 
consider corresponding noninteracting problems, which 
lead to the same four-current and four-potential as the 
original interacting system. Due to the quantum elec- 
trodynamical Rungc-Gross theorem, the system is fully 
described by the four-current and four-potential alone 
and one can calculate in principle all observables by only 
knowing these two quantities. The effective external four- 
potential and the effective four-currents of the Kohn- 
Sham system are uniquely defined. Especially, in con- 
trast to the considerations in the effective external 
current is shown to be simply the sum of the external 
current of the interacting problem and the four-current 
of the matter system. This provides the foundations to 
perform ab initio QED-calculations by solving a corre- 
sponding noninteracting quantum system coupling to an 
effective external potential and an effective external cur- 
rent. For special forms of the initial noninteracting con- 
figuration the noninteracting problem separates into sim- 
ple Dirac-Kohn-Sham equations and a Maxwell equation. 
Finally, assuming the initial configuration to be of this 
special type, we gave a straightforward approximation for 
the effective four-potential of the noninteracting matter 
system. 
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